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We analyze theoretically the effect of friction on quartz 
crystal microbalance (QCM) measurements that probe soft 
(viscoelastic) films and biomolecular layers adsorbed from 
aqueous solutions. While water provides a natural environ- 
ment for biomolecules, an interface with unknown rheological 
properties forms between the adsorbed soft molecular layer 
and the quartz substrate in the latter case. We investigate 
therefore the dynamics of soft films adsorbed onto a solid 
quartz surface within a continuum mechanics approach using 
both the Maxwell and the Voight /Kelvin models of viscoelas- 
ticity and their combination. The rigorous expressions de- 
rived for the acoustic response of a quartz crystal oscillator, 
accounting for both interfacial (sliding) friction and internal 
friction (viscosity) , demonstrate that the QCM can be used as 
a sensor for quantitative characteristization of friction effects 
as well as for "in situ" measurements of mechanical properties 
of adsorbed biomolecular films. 



I. INTRODUCTION 

A universal slip law for a moving interface between 
a solid and a fluid experiencing shear stress, has been 
demostrated recentlylil for cases when the no-slip condi- 
tion breaks down. Experimentally, this situation can be 
addressed with a piezoelectric (quartz) plate, which os- 
cillates in thickness-shear mode while in contact with a 
fluid mediurrn. Traditionally, such quartz crystal oscil- 
lators have been utilized as a microbalance (QCM) for 
weighing a negligible amount of_adsorbed mass in vac- 
uum or in a gaseous environments. Current applications 
of the QCM in a vacuum or gas phase range from nar 
otribology experiments on thin films of solidified gasesl 
to sensor measurements of viscoelastic properties of thin 
films deposited onto the quartz-crystal surfacell In re- 
cent liquid phase experimentsBEB, the QCM technique 
has been adopted as a biosensor for measuring the prop- 
erties of adsorbed polymer- and biomolecular films. How- 
ever, a finite amount of slippage arising from a weak cou- 
pling between film and oscillating substrate can change 
both the resonance frequency / and the quality factor 
Q of the oscillator and renormalize the results of mea- 
surements. Therefore, understanding the role of sliding 
friction in liquid QCM experiments is of key importance. 
An extensive review of the theory of sliding friction and 
the quartz crystal microbalance method can be found in 
the recent book by Perssontil. 

For a quantitive characterization of slippage effects it 
is convenient to use the ratio between the shift A/ of the 
QCM resonance frequency and the inverse quality fac- 



tor A(Q 1 ). This measure was. introduced in the QCM 
literature by Krim and WidomQ as a "slip time" r, 

r = ^2. (i) 

In the following we prefer to refer to the dissipation factor 
D = Q^ 1 rather than to the inverse quality factor of 
the oscillatom. The shifts AD and A/ can be measured 
simultaneously in QCM experiments; one always finds 
AD to be positive while the resonance frequency shift 
A/ is negative since mass loading lowers the oscillator 
resonance frequency compared to the unloaded state in 
vacuum. _ 

According to the results of Krim's groupEl, a partial 
decoupling of the overlayer should occur when 2-7t/t > 
0.5. Typically, a quartz crystal oscillates with frequency 
/ ~ 10 MHz and therefore decoupling starts when r > 
t c ~ 10~ 8 s. Ear instance, in QCM experiments with 
thin water filmsB A/ - 1 - 10 Hz, AD - 10~ 7 - 10~ 6 
and hence r ~ 10~ 9 — 10~ 7 s, which is quite close to t c . 
It follows that a relatively low decoupling threshold can 
lead to noticeable interfacial friction effects even in the 
case of molecular ly thin films. 

In the nanotribology QCM-experiments of Krim and 
collaborators, where they studicjdiiinas of solidified gases, 
condensed onto a solid substrateaEfQu, the interfacial fric- 
tion coefficient 7/ for a thin rigid film was found to be 
inversely proportional to the slip time and proportional 
to the surface mass density to/ = pfhf of the film, 



7/ 



V l T - 



(2) 



This result can also be viewed as the definition of the slip 
time r, which isa the characteristic time it takes for the 
film velocity to decrease by a factor 1/e. However, for soft 
interfaces Eq. (0) becomes invalid because of additional 
viscous dissipation of energy in the material. 

In this article we analyze theoretically the QCM re- 
sponse, accounting for both viscosity and interfacial fric- 
tion effects, and derive how the sliding friction coefficient 
is related to the values of A/ and AD, which can be mea- 
sured. Throughout the paper we use the term "soft" to 
characterize properties which are opposite to "rigid" , i.e., 
the viscoelasticity of the material. 

We will show that under certain conditions, the dis- 
sipation factor AD can have a maximum as a function 
of the viscosity as well as of the sliding friction of the 
QCM-probed material. We derive the exact expressions 
for the maximal dissipation shift. The general analysis 
presented here is valid for arbitrary layer thickness. Slid- 
ing friction of ultrathin films, adsorbed from a gaseous 
(or liquid) phase, emerges as a special limiting case. 



From a practical point of view, one of the particularly 
important topics treated here is concerned with liquid- 
and biosensor QCM applications. Here the oscillator, 
with its adsorbed biomolecular layer, operates in a bulk 
viscous liquid (water), which provides a natural environ- 
ment for in situ measurements. We have modeled this 
system as a layered soft (viscoelastic) material covered 
by a bulk newtonian liquid. Linear viscoelasticity of dif- 
ferent biological interfaces has been analysed for both 
Maxwell and Kelvin/Voight schemes and their combina- 
tion. 



and Stockbridge's result for the QCM dissipation shiftl 



a 




FIG. 1. Schematic depiction of a quartz crystal oscilla- 
tor (microbalance) loaded on one side by a layered medium. 
Layer I could, e.g., be a thin biomolecular layer, layer II could 
be, e.g., bulk water. 



m q V ^ 

Here m q = p q h q is the surface mass density of the quartz 
plate, p q and h q is the mass density and thickness of the 
quartz, respectively. In the general case of a layered vis- 
cous (viscoelastic) medium, both A/ and AD are com- 
plex functions of thickness, mass and viscoelastic prop- 
erties of the materialliirlHl. These characteristics can be 
derived from the acoustic impedance Z(lo) of the system. 
The acoustic impedance 



Z = Z(lu) = Z — iZ 



(3) 



is equal to the mechanical impedance per unit area. The 
real part Z of Z(lu) is a resistive term, which is propor- 
tional to the average energy dissipation per unit time, 
while the imaginary part Z is a reactive term associ- 
ated with the inertia of the oscillatorQ'Q. For a material 
rigidly attached to the quartz surface, the so-called no- 
slip boundary conditions are valid, i.e. the velocity of 
the quartz surface is equal to the velocity of the adjacent 
layer. Below we indicate acoustic characteristics deter- 
mined in the absence of slippage (no-slip boundary con- 
ditions) by a tilde sign (~) over the appropriate symbol. 
The shifts in the QCM resonance frequency and in the 
dissipation factor in the absence of slippage are given by 
the expressions^ 



II. DAMPING OF A QUARTZ OSCILLATOR BY 
SLIPPING SOFT MATERIAL. 

Let us consider a quartz crystal microbalance (QCM), 
oscillating in thickness-shear mode and on one side cov- 
ered by a soft adsorbed layer (as depicted schematically 
in Fig. 1). We describe the motion of the QCM as the 
forceidpvibration of a damped oscillator driven by periodic 
forceffll. 

When mass is deposited onto the crystal surface, the 
resonance frequency changes. For a small amount of mass 
or a molecularly thin film rigidly attached to the QCM 
surface, the shift of resonance frequency A/ is propor- 
tional to the density p and thickness h of the overlayerEJ, 
while there is no dissipation shift AD. When the oscilla- 
tor operates in a bulk liquid environment, both A/ and 
AD are finite and depend on the density p and the vis- 
cosity r\ of the liquid. Classic results are Kanazawa and 
Gordon's resonaace frequency shift of a QCM immersed 
in a bulk liquidlia: 
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and 



AD = Rc 
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(4) 



(5) 



Let us now consider the effect of potential slip arising 
from a weak coupling of film and oscillating substrate. In 
this case the no-slip boundary condition breaks down and 
a finite difference between the velocity qo of the quartz 
surface and the velocity vq of the adjacent layer at the in- 
terface (index 0) appears leading to an interfacial friction 
force, F,, where 



F s =7s(gb -Vq). 



(0) 



Here 7 S is the proportionality coefficient characterizing 
sliding friction. It is a constant in the low shear rate 
approximation. Slippage can also be characterized by a 
slippage coefficient A which is simply the inverse of the 
sliding friction coefficient!!!!, 



A = 7 S 



(7) 



The no-slip condition corresponds to A — » 0, while the op- 
posite limit A — > co gives infinite slippage and a complete 



decoupling of the overlayer from the substrate. Effects of 
slippage can easily be included in the equation of motion 
for the-pscillator. Following the damped oscillator model 
of Ref J£l, we may express the shift of resonance frequency 
and dissipation factor in the presence of slippage as 



A/ 



-Ira 



Z 



2irm q {\ + \Z) 



AD &Re 



7rfm q (l + XZ) 



(8) 



(9) 



where as already mentioned A = 1/7 S . 

Following the approach of Krim's groupOTJ, we intro- 
duce the characteristic slip time r s as the ratio between 
the shift of the dissipation factor and the change in res- 
onance frequency change, 



AD 

AirAf 



(10) 



Proceeding by means of an analogy, we define another 
ratio with dimension of time, 



TO 



AD 

4ttA/' 



(11) 



in order to characterize the temporal QCM response 
when no-slip boundary conditions apply. Defined in this 
way, the characteristic time To can be attributed solely to 
the effects of internal friction (viscosity) on the damping 
of the QCM, i.e, AD is proportional to the resistive term 
(||) in the acoustic "no-slip" impedance Z. 

Using Eqs. (p|)-(p"l|), we may express the slippage coef- 
ficient A and the interfacial friction coefficient 7 S as func- 
tions of t s and tq. One finds 



A 
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and 



j s = 2nAfm q 
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(12) 



(13) 



Let us consider a viscoelastic medium of finite thick- 
ness h in contact with the quartz oscillator (Fig. 1). ft 
is known, that the acoustic response of the oscillator 
strongly depends on the overlayer thickness. "Thin" or 
"thick" films refers to the film thickness being smaller or 
greater than the inverse values of the decay constant a 
and propagation constant k of the-acoustic waves propa- 
gating in the coated quartz plateE3. 

In the limiting case of an ultrathin film (ha -C 1, hk <C 
1), acoustic shear waves can propagate through the ad- 
sorbed film without dissipation. Such a layer has a Sauer- 
brey solid-like response^, for which 

A/=-/o^, AD = 0, 

m„ 



and where mj = phf is the mass per unit area of the 
film. Substitution of A/ and AD in Eq. ([13]) gives the 
coefficient of interfacial friction. For a molecularly thin 
overlayer one finds 

7 s0 = m//r s , 

which reproduces the result of Krim and WidomEl. For 
slipping soft material of arbitrary thickness, one must 
take into account both interfacial friction and viscosity 
and use Eqs. (§)-([l3|). 

An analysis of the general expressions (Q), (||) and (||) 
allows the dissipation factor AD (A) to be expressed as a 
function of slippage as 



AD(A) = 
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(14) 



From Eq. (flj) it follows that the dissipation shift is a 
nonmonotonic function of slippage A, which peaks for a 
critical value A*. If we take the limit Z' = correspond- 
ing to a rigid thin film, the maximal dissipation is 



AD(A* 



f 



(15) 



The dissipation is maximal in this case (wr = 0) when 
the condition 



UJT S 



1 



is satisfied, a criterion for_sliding thin rigid films first 
found by Krim and WidomQ. 

The condition (fi^) is readily understandable. For a 
non-dissipative medium, the damping of the oscillator 
is proportional to the inertial contribution A/// ~ rrif 
associated with the slippage of the added mass pushed 
along the quartz surface. A rigid thin film will give rise 
to maximal dissipation when the sliding friction reaches 
the value 

7* = rrif ■ u> , (16) 

or equivalently when the slippage coefficient is 

X* = l/(m f -cj). (17) 

This maximum value of the dissipation factor is 



AD(A*) = m f /m q 



(18) 



The lower the surface mass density is, the larger is the 
slippage coefficient A* for which the dissipation factor is 
maximal, and the smaller is this maximal value (Fig. 2). 
An important consequence of the result (|8|) is that the 
maximum value of the dissipation factor is frequency in- 
dependent. Therefore, one may consider the presence 
of such a frequency independent maximum as a test of 
whether slippage occurs in dynamic QCM measurements; 



if the dissipation peak remains constant when the oscil- 
lation frequency u> is varied, the peak can be attributed 
to sliding friction. 
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FIG. 2. Sketch of typical nonmonotonic dependence of the 
dissipation factor AD on the dimensionless slippage coeffi- 
cient \ujrrif for a thin, rigid film adsorbed onto the quartz 
plate. The calculation was done for film thickness hf = lfim 
and film mass density pf — 1 g/cm 3 . 

Experimentally, the sliding friction is an external pa- 
rameter which can be changed, e.g., by varying pressure 
or temperature. Maxima in the dissipation has been ob- 
served experimentallyura for solid molecularly thin films 
sliding along a quartz surface oscillating in gaseous envi- 
ronments. The maxima appeared when the the condition 
ljt s = 1 was met. In these experiments, the effect of in- 
terfacial friction on the QCM damping was obvious from 
the fact that viscous losses are negligible both in a gas 
and in rigid molecularly thin films. However, in exper- 
iments with soft films, the dissipation factor can reach 
a maximum as temperature or pressure is varied, since 
these affect the internal friction (shear viscosity) of the 
material. When slippage is absent (or negligibly small), 
viscoelastic effects dominate and the dissipation maxi- 
mum becomes frequency dependent (see Section 3C). 

In case of sliding thin viscoelastic over layers, internal 
friction (viscosity) leads to additional but small contri- 
butions to the QCM damping since AD ^ 0. As a result 
the position of the dissipation maximum is only slightly 
shifted. On the contrary, the viscosity of soft materials 
of finite thickness will significantly influence the damping 
and produce a "viscous" -type dissipation peak. In order 
to separate these effects, we study below the no-slip lim- 
iting case — when the interfacial slippage vanishes — for 
the viscoelastic film dynamics . 



III. ACOUSTIC IMPEDANCE OF A QCM 
LOADED BY A VISCOELASTIC LAYER 

Viscoelasticity is a common feature of complex fluids 
and in polymer rheology. A limiting case is that of highly 
viscous fluids - amorphous solids and glasses. Polymers 
exhibit linear viscoelasticiteiif-both deformations and de- 
formation rates are smal£3'E2l. In the linear viscoelas- 
tic regime, there is proportionality between stored en- 
ergy aed strain and between strain rate and dissipated 
energyEl The most well known mathematical formula- 
tions of linear viscoelasticity are the Maxwell model of 
relaxation in highly viscous fluids and the Kelvin/Voight 
model of viscoelastic solids. 



A. Highly viscous Maxwell fluids 

Let us derive the acoustic impedance of a QCM cov- 
ered by an overlayer of highly viscous material. A very 
viscous complex fluid responds initially to deformation 
in the same way as an elastic solid, but after the defor- 
mation stops the fluid relaxes. This relaxation processes 
can be characterized by the viscous relaxation time tm, 
which is a measure o-Lthe time it takes for the remaining 
stress to be dampedEJ. Such viscoelastic behavior of the 
complex fluid can be characterized by the shear viscos- 
ity coefficient r\ and the shear elasticity modulus /i. The 
characteristic relaxation time tm is of the order of their 
ratiotL2l: 

TM ~ V/fJ- 

Viscoelastic properties of a complex fluid can be 
treated within the Maxwell model which represents a me- 
chanical analogue of viscoelasticity via a dashpot (vis- 
cous newtonian element) and a spring (elastic hookian 
element), arranged in series (Fig. 3a). A Maxwell fluid 
subject to quartz crystal oscillations of frequency lo be- 
haves as a newtonian viscous fluid with shear viscosity 
r\ if tm « w _1 . In this limit, the relation between the 
stress tensor an- and the strain tensor Uik is the same as 
for the ordinary newtonian fluid with a^ = 1uor\Uik- In 
the opposite case of large u> ^$> r _1 , the maxwellian fluid 
responds to shear stress as a solid. The strain-stress reis*- 
tion <jjfe = 2{j,Uik describes the elasticity of a solid bodylia. 



,1, 




FIG. 3. Basic models for the linear viscoelasticity of differ- 
ent materials, a) Maxwell fluid; b) Voight/Kelvin viscoelastic 
solid; c) and d) more complex rheological models of viscoelas- 
tic solids. 



The equation of motion for a maxwellian fluid, which 
includes both the limits of slow- and fast motion, is given 
by the expression 



2fi- 



du 



ik 



(11 



d(Jjk _j_ <Hk_ 

dt 



TM 



(19) 



Equation (R9h assumes that the internal stress is expo- 
nentially damped on the characteristic time scale tm af- 
ter the applied stress is removecO Using Eq. (|19|) we 
find the stress-strain relation for maxwellian fluid to be 
of the form 



Oik 



2pu tk 



1 - i/ujT M ' 
which can be rewritten as 

a ik = 2p*u ik . 



(20) 



(21) 



Here p,* M = p M + ip" m is the complex shear modulus of 
a Maxwell fluid. The real part (storage modulus) is 



' r] z to • p 

Mm = 2 I ? 2 ' 

p Z + U) Zr q Z 

while the imaginary part (loss modulus) is 
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(22) 



(23) 



Using the strain-stress relation (pfl) and the conserva- 
tion of momentum density, 

du ik d 

--d^ = di {pvi) > 



we get a wave equation, 
: d 2 u x (y,t) 



M 



dy 2 



-puj 2 u x (y,t), 



(24) 



for bulk shear waves propagating in the system. Here 
u x (y,t) is the a component of the displacement vec- 
tor (Fig. 1). Solving eauation (|24|) with no-slip bound- 
ary conditions, we findtil the acoustic impedance of the 
maxwellian viscoelastic fluid overlayer of thickness h to 
be given by the expression 



£m*i 



2£.h 



iuj 1 + e 2( - h 



(25) 



Here £ = a + ik] a is the decay constant and k is the 
propagation constant. These quantities are given by for- 
mulae 



OLM 



\\J\I 1 + Xm-Xm (26) 



In Eqs. (|26|) and (|2jU,the viscous wave penetration depth, 
5 = (2rj/puj) 1 ^ 2 isEj the distance over which the am- 
plitude falls off by a factor of e, the viscoelastic ratio, 
X = M 1 1 1 i i s the ratio between real part (storage mod- 
ulus) and the imaginary part (loss modulus) of the com- 
plex shear modulus. For Maxwell material 

xm = Wm- 

Below we will characterize the material by the viscoelas- 
tic ratio x rather then by the relaxation time tm- These 
quantities are related, since xm ~ utm- 



B. Thin Maxwell-fluid over layers 

Here we derive the shifts A/ and AD characteristic of 
a Maxwell fluid thin film rigidly attached to the surface 
of a QCM oscillating in vacuum. By a series expansion in 
ha <C 1 and hk <S^ 1, which is permissible in the special 
case of a thin overlayer, we find from (|2^) - ( P7|) and from 
(|) and (§) that 
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AD « (2h 3 p 2 Lo)/(3m q r]) 



(28) 

(29) 



It should be noted, that Eq. ( |29| ) has the same form as 
the dissipation factor for a thin pure viscous layer with 
constant shear viscosity coefficient r\. 

From expression (|28|) we obtain for the equivalent mass 
rriM of the film 



mm = rn/(l + h 2 pu> 2 /3/z) 



(30) 



which differs from the actual mass rrif and depends on the 
shear elasticity p of the material. Our results (S)-(pO) 
are in agreement with those of Johannsmann et alx3, who 
calculated and measured by QCM the equivalent mass of 
a viscoelastic Langmuir-Blodgett film. The characteristic 
"viscous" time tq for a thin Maxwell film is, therefore 



r « (h 2 p)/(3r,) . 



(31) 



Because tq is defined as the ratio between the experimen- 
tally measured values of AD and A/, one can calculate 
from Eq. ( j3l| ) the viscosity of a thin Maxwell film if the 
film mass density and thickness are known or measured 
separately (e.g., by neutron scattering or ellipsomctry 
methods) . 

Let us now proceed with a comparison between the re- 
sults for a Maxwell fluid with those for a thin viscoelastic 
solid film. 



C. Solid viscoelastic Voight/Kelvin overlayers 



km = t 



1 + Xm + Xm ■ 



In the Voight/Kelvin scheme (Fig. 3b) for describing 
\^'J visoelastic solid materials, a complex viscoelastic shear 



modulus, jjiy — /j,y + ijjl'vt is used. Its real part is the 
storage modulus 

Py = P, (32) 

while the imaginary part is the loss modulus 



p y — uj-q . 



(33) 



The dynamics of viscoelastic solids can alsp-be charac- 
terized by the retardation time Ty ~ r/y / pyBj or by the 
viscoelastic ratio 
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(34) 



The acoustic impedance of a quartz oscillator in the 
absence of slip is given by expression (Gq) where £y = 
ay + iky and 



(35) 
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(36) 



In the same manner as in the previous section, we find 
from (|2q) taken together with (B5|) and (Bq) that 
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(37) 



(38) 



The characteristic time tq depends on frequency and 
viscoelastic moduli of the material. One finds 



to 



h 2 prjuj 2 



3(p 2 +T] 2 L0 2 ) 



We therefore obtain the equivalent mass of the 
Voight/Kclvin layer as 



my = rrif I 1 



h 2 puj 2 p 
3(p 2 + i] 2 uj 2 ) 



(39) 



which includes a correction of the true film mass density 
to/ due to the viscoelasticity of the material. 

Equations (|3q) and ( p9|) together with to allow us to 
calculate the viscoelastic parameters py and r\y from si- 
multaneously measured values of A/ and AD if it possi- 
ble to control the film thickness h independently. 

From expression ( pSf ) it follows that if during the ex- 
periment the shear elasticity is constant but the viscosity 
varies, the dissipation factor as a function of shear viscos- 
ity rj will have a maximum when r/* = p/uj (or xv = 1, 
respectively) (Fig. 4). And vice versa, in experiments 
with constant viscosity but varying shear elasticity p, 



the resonance frequency shift will have a maximum for 
p* = rjuj (xv = 1) (Fig. 5). The maximum value of the 
dissipation factor 



AD(r?*) = 



h 3 p 2 u; 2 
3m q p 



depends on the vibration frequency u> squared. Since the 
dissipation peak due to sliding friction was found to be 
frequency independent, this strong frequency dependence 
of AD(ij*) can be used as a test to see if slippage is 
absent in QCM "viscoelastic" thin film measurements. 
Note that the "viscous" peak value AD(rf) is weaker by 
a factor h 2 pu> 2 / p than the peak value AD(X*) of Eq. JlS}). 
It should be noted, that the frequency behavior of the 
QCM characteristics in the linear viscoelasticity region, 
viz. A/(w), can also test for the type of material - either 
viscoelastic fluid or viscoelastic solid. This is because 
from our results ( p8| ) and (|3g) for A /a/ and Afy, it fol- 
lows that we expect a straight line if A/(w) is plotted vs. 
lo 2 for the case of a Maxwell fluid and a deviation from a 
straight line for other types of viscoelastic materials. In 
the case of Maxwell fluid, the slope of the straight line 
gives the elastic modulus p of the material. (The vis- 
cosity t\m of the Maxwell fluid film can then be deduced 
directly from the measured values of AD or To if the film 
thickness is known.) 




FIG. 4. Calculated dissipation factor AD of a quartz 
crystal microbalance (QCM) covered with a thin film of a 
Voight/Kelvin viscoelastic solid as a function of its inverse 
viscoelastic ratio x _1 — f\^lp\ hi = 1 /im, u) = 2-k ■ 10 MHz, 
p = 10 10 dyn/cm 2 , p = 1 g/cm 3 . 
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Here /ii,pi is the thickness and density of the first layer 
while A/i 2 , p 2 refer to the second layer, 773 and p% denotes 
the viscosity and density of the water solution, u* ^ ,0:1.2 
and fci, 2 are given by Eqs. (§1), (||), ©, and (§3) for 
Maxwell layers or by Eqs. (32) - (pq) for Voight/Kclvin 
layers, respectively. The acoustic response A/ and AD 
can be calculated using the general results (0) and (0) 
together with (Q) and (J5J) , where Z& has to be substituted 
for Z. 

In the next section we analyse formulae (EH) - (E3) in 
the limiting case of a thin viscoelastic overlayer rigidly 
attached to the quartz plate and loaded by water bulk 
solution on top. 



FIG. 5. Relative resonance frequency shift Sf of a QCM 
due to an adsorbed viscoelastic solid thin film; hi =0.1 /im, 
u = 27r ■ 10MHz, rj = 1 dyn-c/cm 2 , p = 1 g/cm 3 , 

5/ = \Af/Afsauerbrey\, Afsauerbrey — ~hp(jj /2TT17l q . 



IV. LIQUID PHASE QCM APPLICATIONS 

Now, to complete the theory, we derive the acoustic 
response of a quartz oscillator covered by a soft bilayer 
when the system is immersed in a bulk water solution. In 
general, the thickness and the mechanical properties of 
the two layers can be different. Such sandwich structures 
are typical for biological QCM experimentalJ'EZfEj where 
the first layer can be a host polymer, an amphiphilic self- 
assembled monolayer or a Langmuir-Blodgett film serv- 
ing as a matrix for the second, tethered polymer- or ad- 
sorbed protein layer or lipid membrane. Water is an im- 
portant component maintaining the native structure of 
biological molecules and biomembranes. 

The no-slip acoustic impedance of a "three-layer" sys- 
tem is given by the formula 



Z h = 



&Mi 1 - /3e 2 ^ lhl 



1-uJ 
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(40) 



where 



P = -7^. ...... S- 1 . . (41) 
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V. APPARENT "DISAPPEARANCE" OF THE 

EQUIVALENT FILM MASS IN LIQUID QCM 

MEASUREMENTS. 

As shown in Section III, the viscosity and elasticity of 
the tested material affect the measured equivalent mass 
m eq of the film, which therefore differs from the "true" 
film mass mj. When the system operates in a solution, 
acoustic shear waves penetrate through the thin overlayer 
and interact with the bulk viscous medium on its top. Let 
us analyse how the viscosity of the solution (index "2" ) 
changes the dissipation and resonance frequency shift of 
the QCM covered by a viscous/viscoelastic film (index 
"1") and how it influences its equivalent mass. 

First, we restrict our attention to the special case of 
a pure viscous overlayer. By a series expansion valid 
for the thin film limit h/8 C 1, we find that the shift 
in resonance frequency when the QCM oscillates in bulk 
water is a function of the mechanical properties of the 
overlayer, p\ and 771, which enter to linear order in the 
(small) film thickness: 
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At the same time, the overlayer contribution to the dis- 
sipation factor is very small and appears only to second 
order in h/S <C 1, 
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where 
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(46) 



is Stockbridge's result for the dissipation of a QCM in a 
purely viscous bulk liquid. We conclude from (H) that 
the equivalent film mass of a purely viscous layer, m ~ 
TO /( 1 — (?72/02)/(?7i/9i), is obviously equal to zero if r?2/02 = 
ViPi- 



For a QCM probing a thin layer of a Maxwell material 
(ha <C 1, hk <C 1) in a bulk solution, we find that the 
resonance frequency shift and equivalent mass take the 
same form as for a purely viscous film. However, the 
dissipation in a Maxwell fluid includes a term linearly 
dependent on the film thickness, which appears due to 
the finite elasticity pi of the material: 



AD 




— ■ (47) 



Finally, we obtain in a similar way the response of a 
QCM covered by a thin Voight/Kelvin overlayer when 
the system operates in a liquid: 
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When r/uj ^i p, viscoelastic solid and viscoelastic fluid 
films dissipate the same amount of energy. In the oppo- 
site case of ijuo ^> p, a viscoelastic solid film will dissipate 
much less energy than will a Maxwell fluid film. 

The equivalent mass of a viscoelastic solid film 



w ~ "1/(1 



mP2Lo m 



pxinl+^r) 2 ) 1 
equals zero when r] = 770, where 
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mP2 

2/91 
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When the viscosity of the top layer vanishes, 772 = 0, vis- 
coelastic corrections to the equivalent mass appear only 
in higher order corrections in the film thickness (see sec- 
tions 3B and 3C). 

A simple analysis of expressions (Eq) and (ff9j) re- 
veals nonmonotonic behavior of the dissipation factor 
AD(p,ri) and the resonance frequency shift Af(p,r]). 
The dissipation factor will reach a maximum when p* — 
j]lu if the viscosity is constant. In the same way, the res- 
onance frequency has a minimum as a function of viscos- 
ity when rf = p/uj (Fig. 6). This circumstance must be 
taken into account in QCM adsorption experiments if vis- 
coelastic parameters and the layer thickness can vary-si- 
multaneously during adsorption from the liquid phaseO. 

It should be mentioned, that the QCM characteristics 
in the case of a Maxwell fluid are monotonic functions of 
film viscoelasticity, decreasing with p and 77, respectively. 




FIG. 6. Relative resonance frequency shift Sf due to a thin 
viscoelastic solid overlayer as a function of its viscoelastic ra- 
tio x = p/ 1 !^- QCM operates in a bulk aqueous solution; 
Sf = |(A/-A/o)/A/s a uer6rej/|, /o is the resonance frequency 
shift in a bulk liquid without overlayer; A/o = —ri/2irm q S, 

AfSauerbrey — ~ hpu}/2irm q . 



VI. DISCUSSION: SOFT MATTER AND 
BIOSENSOR APPLICATIONS OF QCM 

One of the current challenges is to apply the QCM 
technique to probing biological materials in solution: 
polymers and complex fluids, which are dominated by 
viscous effects. Typically adsorbed from a liquid phase, 
they form soft biological interfaces of multilayer archi- 
tecture at the solid substrate. Recently it has been re- 
ported that mechanical properties of (bio)polymer and 
amphiphilic self- assembled films, adsorbed films of pro- 
teins, enzymes, membrane microspheres (vesicles) and 



■iawstigated by using quartz crys- 
|j25|,|2i| j£ seems crucial that an ad- 



even living cells, 
tal oscillatoraHJli 
equate physical model is chosen for a quantitive inter- 
pretation of such QCM measurements to be possible. In 
addition, biosensor applications of the quartz crystal mi- 
crobalance (QCM) technique are faces with the need to 
control the interfacial friction, when the probed biopoly- 
mer material can slip relative to the quartz surface during 
its oscillatory motion. The slip of the adjacent layer, af- 
fecting the resonance frequency shift A/ and dissipation 
factor AD of the oscillator arises from the weak coupling 
between film and substrate. 

The limit of strong coupling of the layer to the solid 
surface corresponds to the widely used "no-slip" bound- 
ary condition. Experimentally, this boundary condition 
applies to thin solid films evaporated onto the quartz 
surface. Among biomaterials, the no-slip conditions can 
easily be realized for Langmuir-Blodgett films, which can 
be extremely strongly coupled to the solid substrate. 

The usage of QCM for biological and biomedical pur- 
poses brings out the problem of accounting for the com- 
plex rheology of the biological interface formed by the 
adsorbed biomolecular films. In some special cases or 



for particular type of shear deformations, the linear vis- 
coelasticity of polymers can be treated within either 
the Maxwell or the Voight/Kelvin schemes. The results 
of the Voight/Kelvin model can readily be applied to 
quartz crystal acoustic measurements of adsorbed pro- 
teins, which conserve their shape during adsorption and 
do not flow under shear deformation as well as to polymer 
films far from the glass transition region. The Maxwell 
fluid scheme is appropriate for polymer solutions that 
are pure liquids (at low shear rates)Ef] and polymers in 
the amorphous state and in the vicinity of the polymer 
liquid-glass transition. Viscoelastic materials with more 
complex rheology (e.g., cells, membranes, liquid crystal 
polymers etc.) can be described by a combination of 
these two basic viscoelastic schemes. Two such schemes 
are illustrated in Figs. 3c and 3d. They can be analysed 
in a similar as we have done here way by substitution 
of the corresponding expressions for the complex shear 
modulus [A* in the general results ( |25| ) - (J27]) and then in 
@ " ©■ For instance, the complex shear modulus for 
the material described by scheme of Fig. 3d, is given by 
the expression 
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(50) 



Such a model well describes, e.g., the viscoelasticity of 
the lyotropic lamellar phase (the so-called "onion phase" ) 
which organizes itself into multilamellar vesicles which 
are closed-packed and fill up the spaceESl. 

It is convenient to describe the frequency behavior of 
a complex viscous medium exhibiting a spectrum of re- 
laxation times in terms of the three-parameter model 
rj* = 77/(1 + iuT) b of the complex viscosity suggested 
by Reed, Kanazawa and KaufmantJ. In this model, r\ is 
the viscosity of the liquid, 6 is a characteristic exponent 
related to the distribution of relaxation times r. In the 
same manner, we can introduce in our model a complex 
shear modulus of the form 

M* = M/(! - i/ur) b 

which is a generalization for the case of a relaxation time 
distribution characterized by the exponent b. In partic- 
ular, the Maxwell fluid corresponds to the value 6=1. 

In summary, recent results show that in cases where 
the adsorbed material forms a soft interface to the quartz 
crystal microbalance (QCM), both internal and interfa- 
cial friction mechanisms may contribute to energy dissi- 
pation. The combined action of these mechanisms may 
cause the variation of the shift AD in the dissipation 
factor to be a nonmonotonic function of material param- 
eters, which can be varied in an experiment. Depending 
on the system, one or the other type of friction may dom- 
inate. We have suggested here a dynamic "slip/no-slip" 
test that can be applied to QCM measurements with thin 
adsorbed films. The test is based on the frequency de- 
pendence of a peak in the dissipation factor. We have 



shown that the peak value of AD caused by sliding fric- 
tion is frequency independent while the viscoelasic peak 
is strongly dependent on frequency. 

Besides being able to understand the role of interfacial 
friction, an important problem for the application of the 
QCM technique to biosensoring remains. How can the 
device recognize what type of material is being probed 
except by detection of an "anonymous" deposited mass? 
In this article we suggest that analysing the frequency 
behavior of QCM characteristics as one of possible test 
what material is overlayered the quartz surface. 

Another way is a direct biochemical modification of 
QCM surface for selective binding of tested biological 
molecules. Recently, an immunosensing system based on 
antigen-coated QCM has been reported elsewhereEH. It 
has been shown that QCM results were in agreement with 
those of enzyme linked immunosorbent assay (ELISA). 
Perhaps the successful road to QCM usage as a biosen- 
sor device, is to combine both of dynamic measurements 
and immunology treatment directions. 

Presented here rigorous formulae for experimentally 
measured A/ and AD values, make it possible to investi- 
gate viscoelastic contribution in a wide range of biomolec- 
ular materials and can be useful for the correct interpre- 
tation of biologically oriented QCM experiments. 
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